Enhanced Pairing Correlations near Oxygen Dopants in Cuprate Superconductors 
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Recent experiments on Bi-based cuprate superconductors have revealed an unexpected enhance- 
ment of the pairing correlations near the interstitial oxygen dopant ions. Here we propose a possible 
mechanism - based on local screening effects - by which the oxygen dopants do modify the electronic 
parameters within the Cu02 planes and strongly increase the superexchange coupling J. This en- 
hances the spin pairing effects locally and may explain the observed spatial variations of the density 
of states and the pairing gap. 
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In terms of global phase behavior, the cuprate fami- 
lies hosting high-T c superconductivity (SC) are similar: 
parent compounds are insulating antiferromagnets with 
Tjv ~ 300 — 400iv , and become SC once doped charge 
carriers destroy magnetic order and fermionic bands are 
developed. This universality manifests that the inter- 
actions driving magnetism and superconductivity are 
mostly confined to the Cu02 planes common to all the 
high-T c cuprates. What is remarkable, however, is that 
the SC transition temperature T c broadly varies from one 
family to another, ranging from ~ AQK in L^-zSr^CuO^ 
to the values as high as ~ 140K in Hg-based compounds. 
A number of ideas concerning the material dependence 
of T c have been discussed: the influence of apical oxygens 
1 1 , charge ordering effects [2J , band structure variations 
3j, the type of dopant induced disorder [4[, etc. Given 
the high sensitivity of correlated electrons in oxides to 
external perturbations, all these material specific factors 
may indeed have a strong impact on how the competition 
between SC and magnetism is resolved in a particular 
compound and thus on the T c values. 

Local disorder - an inevitable side effect of chemical 
doping - triggers also nanoscale phase inhomogeneities 
commonly observed in cuprates. Signatures of the pair- 
ing effects well above the bulk T c (see, e.g., 0-0]) nave led 
to a notion of " local pairing gaps" and " local T c " . Typ- 
ically, the dopant ions are expected to perturb the lat- 
tice and chemical bonds, reduce the electronic mobility, 
and thus suppress SC in favor of magnetism Q ■ Strong 
enhancement of the pairing correlations near the inter- 
stitial oxygen dopant ions in Bi 2 Sr 2 CaCu 2 8+ ,5 reported 
by McElroy et al. Q was therefore a big surprise. Their 
scanning tunneling microscopy (STM) data have been 
nicely reproduced by the model calculations @ assum- 
ing a positive impact of the dopant ions on the pairing 
potential, and triggered a broad discussion 1(J- H| on the 
origin of this striking observation. 

Given the essential role of magnetic interactions in 
cuprates, there has been a strong focus on the possible 



increase of the spin exchange coupling J by dopant ions. 
In principle, this "quantum-chemistry" parameter is sen- 
sitive to local electronic structure and widely varies in 
cuprates, from ~ 110 — 140 meV in SC cuprates to the 
values as high as ~ 240 meV in a quasi-one dimensional 
compound Sr 2 Cu03 (seeRef. [16| for related discussions). 
However, the most detailed calculations 3, 14 1 found so 
far only moderate impact of dopants on J. 

In this Letter, we propose a mechanism by which the 
dopant oxygens may indeed strongly enhance the inter- 
action J locally. Our key idea is to go beyond static level 
structure used in standard calculations of J, by taking 
into account a dynamical change of electronic parameters 
due to polarization effects. Different from previous work 
considering a dopant ion as mere point charge modifying 
the Madelung potentials on Cu0 2 planes (related changes 
in J are small l_3j, llj] ) , we notice that there is a strong 
covalency between the dopant Od and closely located api- 
cal Q oxygen electrons, forming a molecular orbital com- 
plex (see Fig. 1). In fact, the mixing of and apical O a 
orbitals is suggested by experiment fl7| . Consider now 
virtual p-d and d-d charge transitions within the Cu-O- 
Cu bond that lead to the spin exchange J. We will show 
that the corresponding excitation energies A pc i and U 
are dynamically screened by the polarization of molec- 
ular orbitals hence enhancing J. The ionic polarization 

here, the 



effects on excitation energies are known [18W2C 
effect is greatly amplified due to cooperative response of 
the spatially extended O a — O^— O a complex having much 
higher polarizability than that of constituent ions alone. 
We will also show, by an exact diagonalization of the 
t — J model, that local enhancement of J leads to the 
spatial variations in density of electronic states (DOS) 
observed in STM experiments 0]- Our findings suggest 
an interesting possibility of quantum-chemistry control 
of the key interaction J in cuprates. 

To begin with, we recall that a dopant oxygen O^ is 
located between the SrO and BiO planes [3, Eo] , at the 
distance ~ 2.2 A from two apical O a ions. In the ground 



2 




(b) 



ft'o,, Pi 



n -h- 

0„ O a2 ft' 



E-E„ 



(C) 



Cu, 



.. +e 




+e 



— • 

Cu 2 Cu 



o 



Cu, 



FIG. 1: (color online), (a) The O a — Od — O a molecular com- 
plex formed by a dopant Od and its two apical O a oxygen 
neighbors. 2p and 3s orbitals involved in the polarization pro- 
cess are depicted, (b) The level structure and hopping path- 
ways of these orbitals. (c) Charge disproportionation due to 
the p-d (left) and d-d (right) electron transfer within the Cu02 
plane. The corresponding excitation energies A p d and U are 
screened by the high-energy transitions between O a — Od — O a 
molecular orbitals. 



state, i.e., before the virtual charge excitation within the 
Cu02 plane is made, Hamiltonian of the O a — Od— O a 
molecule comprises two terms, H mo i 
where the first term 



H ion ~\~ H COVl 



Hion = [E { p d) n p + E^n.] + £ [E p n^ + E s n^} (1) 



stands for the 2p and 3s electron energies of a dopant 
(E p d }) and apical (E PtS ) oxygen ions: n p — p^p and 
n s = s's are the corresponding particle numbers. Ex- 
perimental data [H and band structure calculations 
(lOJ suggest that 2p level of a dopant is higher than 
that of apical O a by several eV's; we will take below 
E p — E p — 2 eV (the results are not very sensitive to 
this parameter). In Eq. ([1]), summation over the spin di- 
rection and 2p orbitals lying in the xy plane formed by 
the O a — Od— O a molecule, see Fig. 1(a), is implied. 

The second term in H mo i describes the formation of 
the molecular orbitals due to the hopping between the 



dopant and apical oxygen states: 

H C ov = - T ppa [p x (upi -vp 2 ) + p\{up 2 - vpi)} (2) 

- T pp7I [pl(up' 2 + vp[) + pj,(upi + vp' 2 )] 

- T sp [pl(usi - vs 2 ) + p\{us 2 -vsi)} 

- T sp [s\ Pl +p 2 )} - T ss [s\ Sl + s 2 )} + h.c, 

where u = cos /3, v = sin/3, and /3 is the angle between the 
apical pi and dopant p x orbitals [see Fig. 1(a)]. Strong 
mixing of electronic states via H cov enhances the polariz- 
ability of the O a — O^— O a molecular complex, the effect 
which goes well beyond a conventional dipolar polariza- 
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tion of oxygen ions alone 
We use the relation T, 



2 - L ppa 



note Tpp a = Tpp. The signs in Eq. @ are dictated 
by the orbital structure in Fig. 1(a), with a convention 
T pp ,T sp ,T ss > 0. Since 3s orbital is more extended than 
2p one, a relation T pp < T sp < T ss holds [HJ. We as- 



L pp 

sume the scaling (T pp ,T sp ,T ss ) = 



T pp (l,K, k 2 ) with a 



representative value k = 4/3, and vary T pp within the 
1.5 — 2.0 eV range. This should be a reasonable esti- 
mate for the overlap between rather extended 2p oxygen 
orbitals at distance R(O a — Od) — 2.2 A; for compari- 
son, the hopping t p d a between the less extended copper 
d 3x 2__ r 2 and in-plane oxygen 2p x orbitals at similar dis- 
1.9 A is about 



22|. 



tance R(Cu - O) ~ 1.9 A is about -1.8 eV 

The exchange J between Cu spins is realized via the 
virtual hoppings of electrons in the Cu02 planes. The 
charges dynamically generated on Cu and O sites during 
these transitions [see Fig. 1(c)] modify the energy levels 
of 2p and 3s electrons on the apical and dopant oxygens: 
4?i(¥>) = E p , s + ip(r a ) and E$ (tp) = E p d) s + y(f d ), 
correspondingly. The energy shifts y(r) are determined 
by Coulomb potentials of the virtual charges. For in- 
stance, the p-d hopping process creates an extra electron 
(hole) on Cu (O) sites [Fig. 1(c) left]. The potentials 
ip(f a ),ip(f d ) of these charges cause a dynamical recon- 
struction of the molecular orbitals and their populations, 
i.e., polarize the O a — Od— O a complex, which results in 
the renormalization of the virtual excitation energies. 

Under the potential ip(r) of virtual charges, the Hamil- 
tonian of the O a — Od— O a complex becomes H mo i(ip) — 
Hion(<p) + H cov , where H ion (ip) = H ion + H v with 

H v = <p(f d )(n p + n s ) + ]T ^(ra)(n p a) +n<°>) 



0=1,2 

M V( P(ra)[si(Pat+P a j) + h: 

a=l,2 



(3) 



The first two terms stand for the energy shifts of 2p and 
3s levels, while the gradient term describes the ionic po- 
larization of apical oxygens (which are close to the Cu02 
plane). Here, fi — (2p\r\3s) is the dipolar matrix ele- 
ment and the unit vectors i and j select a proper 2p or- 
bital oriented along the gradient V</?. The value of \i can 
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be estimated from the polarizability a of O 2 
0.5 



ion [18j|: 

0.6 A, considering a ~ 2 A 3 
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±y/a(E s - E p ) 
23] and the 3s-2p level separation B s — i5 p ~ 7 — 10 eV 
24j . The gradient term has been included because it 
may change the orbital populations; however, this purely 
ionic effect is found to be small. 

We assume that the electronic transitions within the 
O a — Od— O a complex are faster than charge fluctuations 
in Cu-O-Cu exchange bonds. Indeed, the (inverse) time 
scale for the polarization process is determined by 3s- 
2p level separation ~ 10 eV. This is much larger than 
the Cu-0 charge-transfer energy A p d ~ 3 eV which is 
the dominant parameter dicta ting the strength of Cu-Cu 
virtual hopping t = t 2 v J A p d [25J and thus J. Approxi- 
mation of an " instantaneous" polarization simplifies the 
calculations enabling us to express the screening effects in 
terms of the energies of O a — O^— O a complex in the initial 
[<p(r) = 0] and intermediate [tp{r) ^ 0] states, i.e., before 
and after the electron hopping in the Cu02 plane is made. 
For a given transition, e.g., Cu-0 p-d one, we consider the 
corresponding Coulomb potentials <p(r) on 0^2 and 
ions and calculate the energies E^ ol — (H mo i(<p) — H mo i) 
and E\° l n = (Hi on (ip) — Hi on ). These quantities represent 
the polarization corrections to the p-d charge- transfer 
energy, with the O a — 0<j— O a covalency effects included 
or only considering the ionic polarization, correspond- 
ingly. The difference 5E = E^° ol — Ef°^ is precisely what 
we are looking for: the renormalization of A p d due to 
high-energy response of the O a — 0^— O a complex. We 
find E^ ol Ef° n ', i.e., the effect is indeed dominated by 
the polarization of molecular orbitals. 

The coupling J is estimated from a familiar expression: 



J = 
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Three distinct intermediate states, generated by p-d (A p d 
process), d-d (U process), and so-called pair hopping 
(2A p( i process with two holes on oxygen) transitions, are 
involved here. For each state, we use the corresponding 
potential tp(f) and evaluate the polarization corrections 
to the excitation energies. 

As shown in Figs. 2(a) and 3(a), the dopant has a sur- 
prisingly strong impact on A p d and U values (however, 
the correction to the 2 A p d— process is negligible). The 
effect is sensitive to the hopping T between the apical 
and dopant orbitals and to the energy E s — E p between 
an empty and occupied states: both parameters directly 
control the strength of covalency within the O a — O^— O a 
molecule and thus its polarizability. It is interesting to 
note that doping by the larger size ions like sulfur S 
should increase the covalency and enhance the effect. 

Having obtained the polarization corrections to excita- 
tion energies, we may address the J value, Eq. (dJ). We as- 
sume that p-d hopping t p d is not affected by dopants [26| , 
and calculate the ratio J / Jq of the spin coupling J on the 
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FIG. 2: (color online), (a) Reduction of the A pt j and U excita- 
tion energies due to a polarization of the O a — Od— O a molec- 
ular complex, as a function of hopping T pp . (b) Enhancement 
of a virtual hopping t = tp d /A p d and the superexchange cou- 
pling J on the Cu-Cu bond below a dopant oxygen. 
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FIG. 3: (color online), (a) Reduction of the A pt j and U exci- 
tation energies as a function of 3s- 2p energy level separation 
E s — E p . (b) Enhancement of t and J values on the Cu-Cu 
bond below a dopant ion. 
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FIG. 4: (color online). Exact diagonalization results for the 
local DOS (per spin direction) on different sites i (=0,1,...) of 
the t — t' — J model cluster shown in the inset. The periodic 
boundary conditions are applied. The next-nearest-neighbor 
hopping t' = — OAt. One particular exchange bond indicated 
by an arrow has an enhanced J = 1.7 Jo (compared to Jo = 
OAt in the "bulk"). On the corresponding sites i — or 
i — 15, the spectral weight is suppressed at the Fermi level 
and transferred to the higher binding energies. The result for 
J — Jo case is represented by a dotted line. 



bond below a dopant oxygen to bare Jo in the "bulk" , 
using A p d = 3 eV, U = 8 eV, and U p = 4 eV as the un- 
screened values. The results are presented in Figs. 2(b) 
and 3(b), together with t/t — A p d/{A p d + SA p d) val- 
ues (note that SA p d < 0). A dramatic, by a factor of 
1.5 or even larger, enhancement of the interaction J is 
observed, mainly due to an increased effective hopping 
t oc 1/ Apd (note that J oc t 2 ). 

The above results demonstrate a pronounced effect of 
the oxygen dopants on magnetic interactions which are 
believed to be relevant for the pairing in cuprates. To 
make a closer link with the experiment, we have calcu- 
lated a local DOS by an exact diagonalization of the t — J 
model (on 20-site cluster with 2 holes) where a particu- 
lar bond has an enhanced J value [53] ■ Reduction of the 
DOS at the Fermi level and strong shift of the hump fea- 
ture to higher binding energy is found for sites with large 
J coupling (see Fig. 4). This is consistent with the nor- 
mal state data by Pasupathy et al. Q, and has a clear 
interpretation in terms of an enhanced spin singlet pair- 
ing near the dopants. Based on the model calculations by 
Nunner et al. [9( , it is natural to think that the same cor- 
relations are also responsible for the local enhancement 
of the pairing gaps in the SC state 0, Q; however, this 
remains to be clarified since our small cluster calculations 
cannot directly address the SC state properties. 

To summarize, the oxygen dopants have a strong pos- 
itive impact on magnetic interactions J and enhance the 
pairing effects locally. The effect is of an entirely dy- 
namical nature and thus goes beyond the band structure 



calculations. Apart from the STM experiments, there is 
an interesting possibility to test our theory by the res- 
onant inelastic x-ray scattering: it monitors the high- 
energy spin excitations [HI, and a broad distribution of 
J values predicted here can be directly observed. 
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